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Abstract. We report EVN, MERLIN and VLBA observations at 18 cm, 6 cm and 3.6 cm of the Seyfert galaxies NGC 7674, 
NGC 5506, NGC 2110 and Mrk 1210 to study their structure and proper motions on pc scales and to add some constraints on 
the many possible causes of the radio-quietness of Seyferts. The component configurations in NGC 7674 and NGC 2110 are 
simple, linear structures, whereas the configurations in NGC 5506 and Mrk 1210 have multiple components with no clear axis 
of symmetry. We suggest that NGC 7674 is a low-luminosity compact symmetric object. Comparing the images at different 
epochs, we find a proper motion in NGC 7674 of (0.92 ± 0.07) c between the two central components separated by 282 pc 
and, in NGC 5506, we find a 3 it upper limit of 0.50 c for the components separated by 3.8 pc. Our results confirm and extend 
earlier work showing that the outward motion of radio components in Seyfert galaxies is non-relativistic on pc scales. We 
briefly discuss whether this non-relativistic motion is intrinsic to the jet-formation process or results from deceleration of an 
initially relativistic jet by interaction with the pc or sub-pc scale interstellar medium. We combined our sample with a list 
compiled from the literature of VLBI observations made of Seyfert galaxies, and found that most Seyfert nuclei have at least 
one flat-spectrum component on the VLBI scale, which was not seen in the spectral indices measured at arcsec resolution. We 
found also that the bimodal alignment of pc and kpc radio structures displayed by radio galaxies and quasars is not displayed by 
this sample of Seyferts, which shows a uniform distribution of misalignment between 0° and 90°. The frequent misalignment 
could result from jet precession or from deflection of the jet by interaction with gas in the interstellar medium. 
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1. Introduction 

Active galaxies are conventionally divided into radio-quiet 
and radio-loud categories. Kellermann et al. (1989) found that 
quasars separated into these two classes when the ratio of their 
radio to optical continuum was considered, and a similar result 



has been found by most (e.g. Peacock et al. 1986, Miller et al. 
1990, Ivezic et al. 2002) but not all (White et al. 2000) work- 
ers since. Members of the radio-loud category include BL Lac 
objects, radio galaxies and radio-loud quasars. Radio-loud ob- 
jects have large-scale radio lobes and compact luminous cores 
which often exhibit apparent superluminal motion. Radio-quiet 
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objects include Seyfert galaxies and radio-quiet quasars; most 
of their radio emission is confined to the sub-kpc scale. 

Recent Very Long Baseline Interfere me try (VLBI) has 
shown that many Seyferts also contain compact nuclear ra- 
dio components with brightness temperatures » 10^ K despite 
having radio luminosities over a factor of 10^ less than the most 
powerful active galaxies (e.g. Ulvestad et al. 1998; Ulvestad 
et al. 1999a; Ulvestad et al. 1999b; Roy et al. 2000; Mundell 
et al. 2000), and that they harbour compact components mov- 
ing at non-relativistic speeds. In all cases, the central masses 
are inferred to be black holes of milUons to bilUons of soliir 
masses, so the central potentials are relativistically deep. Since 
most accretion energy is released close to the hole, jets might 
be expected to be relativistic in all cases. Why then are the ra- 
dio sources in Seyferts small, weak and slow compared to those 
in radio-loud objects? 

Various answers to this question have been proposed, in- 
volving either "intrinsic" differences in the central engine or 
"extrinsic" differences in the surrounding medium. Intrinsic 
differences that have been suggested include (i) systematically 
lower black hole masses in Seyferts (Laor 2000), (ii) lower 
black hole spins (Wilson & Colbert 1995), (iii) a "magnetic 
switch" that was identified by Meier et al. (1997) during numer- 
ical modeling of jets, (iv) the production in Seyferts of buoyant 
plasmons that bubble up through the density gradient of the 
narrow Une region (NLR) instead of a collimated relativistic jet 
(Pedlar et al. 1985, Whittle et al. 1986, Taylor et al. 1989), (v) 
a large thermal plasma fraction in the jet in Seyferts (Bicknell 
et al. 1998), or (vi) radiative inefficiency (Falcke & Biermann 
1995). Extrinsic differences generally invoke the rapid decel- 
eration of initially relativistic jets by collisions in a dense sur- 
rounding BLR or ISM (e.g. Norman & Miley 1984). With so 
many possible causes, and so few observational constraints, the 
question remains open why Seyfert galaxies have such low ra- 
dio luminosities. 

One approach to investigate whether the radio quietness of 
Seyferts is due to intrinsic or extrinsic causes is to measure 
speeds of radio components as close to the nucleus as possible, 
hopefully before interactions with the ISM affect them much. 
How close that needs to be is unclear, but signs of jet-NLR 
interaction have been seen on scales of tens of pc (Whittle et al. 
1986; Wilson & Tsvetanov 1994; Bower et al. 1995; Capetti 
et al. 1997; Falcke et al. 1998), so sub-pc scales are desirable. 
Proper motions of Seyfert galaxy radio components have now 
been measured with VLBI in a few cases (e.g. Mrk 231 and 
Mrk 348 by Ulvestad et al. 1999b, NGC 1068 by Roy et al. 
2000, and NGC 4151 by Ulvestad et al. 1998, Ulvestad et al. 
2002), and the motions were found to be < 0.25 c. In III Zw 2, 
V < 0.04 c was seen between two barely-resolved components 
over eight months, followed by a sudden displacement at v > 
1.25 c over seven months (Brunthaler et al. 2000). 

The observations presented here were made for three 
reasons: 1) to measure proper motions of radio components 
on pc scales. 2) to look for extended, low surface-brightness 
radio emission on the scale of the NLR. In the event that slow 
proper motions were measured, the next question would be 
whether components were ejected slowly or whether they 
were ejcted at high, perhaps relativistic velocities and then 



slowed by interaction with the NLR gas. Were such braking 
to occur, one would expect shocks in the jet and possibly jet 
disruption. The in-situ particle acceleration associated with 
shocks could lead then to extended regions of radio emis- 
sion. 3) to measure spectral indices, which provide a useful 
diagnostic of emission mechanisms, source compactness, and 
possible foreground free-free absorption. Most spectral index 
measurements of Seyfert galaxies to date have used the VLA, 
ATCA or WSRT, whose relatively large beams can average 
together many components in the nucleus. Measurements 
with parsec-scale resolution in Seyfert galaxies enable one 
to separate a possible compact, flat-spectrum component 
from surrounding steep-spectrum, optically-thin synchrotron 
emission, but such measurements require VLBI with matched 
beams at two frequencies and are still relatively few. In 
contrast, powerful radio sources have been well studied with 
VLBI and flat-spectrum core components are often seen at the 
base of a conically-expanding jet, thus indicating the location 
of the nucleus. Do Seyfert nuclei also show a single, compact, 
synchrotron self-absorbed core component at the base of a 
coUimatedjet? 

To meet these goals, we observed using Global VLBI, the 
European VLBI Network (EVN), the Multi-Element Radio 
Linked Interferometer Network (MERLIN) and the Very Long 
Baseline Array (VLBA) the Seyfert galaxies NGC 5506 and 
NGC 7674 at multiple epochs for proper motion measurements 
of nuclear radio sources, and Mrk 1210 and NGC 21 10 for the 
first epoch. To enlarge the sample size for deriving statistics on 
the spectral and geometric properties of the pc scale structure of 
Seyferts, we made a complete literature review of VLBI obser- 
vations of Seyfert galaxies. We assume Hq = 75 km s'' Mpc ' 
and qo - 0.5 throughout this paper, and calculate spectral in- 
dices using S ocv°'. 

2. The Sample and Observations 

The galaxies presented here are part of a sample selected for 
global VLBI observations in 1994 to investigate the dichotomy 
between radio-loud and radio-quiet AGN. The sample com- 
prised 16 Seyfert 2 galaxies from Veron-Cetty & Veron (1991) 
with 6 cm core flux densities > 70 mJy and S > -10°. The lit- 
erature was then carefully reviewed for updated spectral clas- 
sifications, some of the objects were discarded and others were 
incorporated based on their new classification and flux density 
measurements. 

Global VLBI, 1994. The sample was observed in 1994 
with global VLBI at 6 cm, using the ten VLBA stations, 
the phased VLA, Effelsberg, Onsala, Medicina and Noto for 
36 min in six snapshots over 1 1 h, yielding a ~ 30 min integra- 
tion of each object (Krichbaum et al. 1994, unpublished). The 
observations yielded only the detection of point sources in the 
cases of NGC 7674, NGC 21 10 and Mrk 1210, and lacked short 
baselines and spectral index information. Only the NGC 5506 
data are included in this paper. 

EVN and MERLIN, 1999/2000. Four of those objects 
(NGC 5506, NGC 7674, NGC 2110 and Mrk 1210) were re- 
observed in 1999 and 2000 with better sensitivity and better 
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Source 


RA (J2000) 


DEC (J2000) 


D 


Type 


Purpose 


Position reference 


NGC 5506 


14 13 14.87926 


-03 12 27.6514 


24.7i«) 


Sin"-) 


Science Target 


Roy et al. (1997)i'^> 


1404-015 


14 04 45.8949 


-01 30 21.937 






Phase calibrator 




NGC 7674 


23 27 56.712 


+08 46 44.14 


116^^' 


Slh2« 


Science Target 


Krichbaum et al. (\99Af'^ 


2327+096 


23 27 33.5808 


+09 40 09.460 






Phase calibrator 


Browne etal. (1998) 


NGC 2110 


05 52 11.376 


-07 27 22.52 


31.13"' 


Sli^w 


Science Target 


Ulvestad & Wilson (1983)3^> 


0541-056 


05 41 38.0848 


-05 41 49.395 






Phase calibrator 




Mrk 1210 


08 04 05.856 


+05 06 49.83 


53.9^' 


Slh"*' 


Science Target 


Falcke et al. (1998)'*"' 


0803+043 


08 03 56.4444 


+04 21 02.724 






Phase calibrator 


Browne et al. (1998) 



Table 1. List of observed sources. Column 2 and 3 give the coordinates used in the observations and for correlation, column 4 gives the source 
distances in Mpc. Column 5 gives the classification according to Veron-Cetty & Veron (2001). Remarks: la) Keel (1996), lb) Nagar et al. 
(2002b), Ic) Coordinates from Roy et al. (1997) observations (unpublished), which based on VLA observations by Ulvestad & Wilson (1984), 
2a) Nishiura et al. (2000), 2b) Veron-Cetty & Veron (2001), 2c) VLA 6 cm (unpublished), 23 25 24.410 +08 30 12.60 (B1950) precessed to 
J2000 by NRAO's scheduling software SCHED, Keel (1996), 3a) Nelson & Whittle (1995), 3b) Veron-Cetty & Veron (1998), 3c) 05 49 46.376 
-07 28 01.99 (B1950) precessed to J2000 by SCHED, 4a) de Vaucouleurs et al. (1991), 4b) Veron-Cetty &. V6ron (1998), 4c) Based on 6 cm 
VLA observations from Falcke et al. (1998). 



(m, v)-plaiie coverage using EVN and MERLIN at 18 cm and 
6 cm. 

The observations were made with eight EVN stations 
(Effelsberg, Cambridge 32 m, Jodrell Bank Mk 1 at 18 cm 
and Jodrell Bank Mk 2 at 6 cm, Medicina, Noto, Onsala 25 m, 
Westerbork array, Torun) at 18 cm (1.65 GHz) on November 
11, 1999, and at 6 cm (4.99 GHz) on February 28, 2000. The 
18 cm EVN observations were recorded in Mklll mode with a 
bandwidth of 56 MHz, 1 bit sampling and dual circular polar- 
ization, and were correlated in Bonn on the Mklll correlator. 
The 6 cm EVN observations were recorded in the new MklV 
mode with a bandwidth of 32 MHz, 2 bit sampling and dual cir- 
cular polarization, and were correlated at the Joint Institute for 
VLBI in Europe (JIVE). The MERLIN observations were cor- 
related in real time with a bandwidth of 32 MHz, 2 bit sampling 
and dual circular polarization. As all sources are weak, phase- 
referencing was used for initial phase calibration. Details of the 
calibrators are given in Table 1 . All calibrators lie within 6° of 
the targets. Phase calibrator scans were obtained every 5 min to 
10 min at 18 cm and every 4 min at 6 cm. 

Additional data. Further observations presented in this pa- 
per were done in the framework of other projects: 

- NGC 5506 was observed during two further VLBI experi- 
ments, designed for proper motion measurement. 

(i) February 5, 1997: NGC 5506 was observed with the 
ten VLBA stations at 18 cm, 6 cm and 3.6 cm with on- 
source times of 91 min, 64 min and 43 min, a bandwidth of 
32 MHz, 2 bit sampling and right circular polarization. 

(ii) May 31, 2000: NGC 5506 was observed with the ten 
VLBA stations at 18 cm and 6 cm with on-source times of 
204 min at each wavelength, a bandwidth of 32 MHz, 2 bit 
sampling and left circular polarization. 

The 1994 and 2000 VLBA observations and the EVN ob- 
servations were reduced in AIPS using the same procedure 
for each (see below), including the same weighting, taper 
and restoring beam size, and the 1997 VLBA observations 
were reduced using the same methods except for slightly 
different tapering. All 18 cm and 6 cm images are shown in 
Fig. 1. 



- NGC 7674 was observed at 18 cm on April 13, 1985 by 
Unger et al. (1988) with four EVN stations using a band- 
width of 28 MHz and dual circular polarization. We re- 
reduced those data from the archive in Bonn following the 
same procedures as for the 1999 EVN 18 cm data, yield- 
ing images with rms - 0.46 mJybeam"' and a beam size 
of 39.3 X 24.2 mas. Amplitude calibration data were not 
archived, so we used T^y^ and gain measurements from the 
immediately preceding experiment (Porcas, priv. comm.). 

- Mrk 1210 was observed on August 6, 1998 with the VLBA 
at 18 and 6 cm with on-source times of 160 min and 
140 min, a bandwidth of 32 MHz, 2 bit sampling and dual 
circular polarization. Phase-referencing was made to the 
same calibrator as for the EVN and MERLIN observations. 
The observations yielded thermal noise limited images 
with resolutions of 18.6 x 14.2 mas and 3.9 x 1.8 mas at 
18 cm and 6 cm. 

Data reduction was done using the Astronomical /mage 
Processing System (AIPS). The EVN and VLBA amplitudes 
were calibrated using Tsys measurements and the MERLIN 
amplitudes were calibrated using the flux density calibrator 
DA 193 assuming flux densities of 2.12 Jy at 18 cm and 
5.15 Jy at 6 cm. Residual delays and phases were determined 
by fringe fitting on the phase reference calibrators and in- 
terpolating the solutions onto the targets. If the target was 
visible in the resulting phase-referenced image, one to three 
iterations of phase self-calibration and deconvolution using 
CLEAN were performed until a thermal noise limited im- 
age was achieved. Absolute position information is then pre- 
served as the self-calibration uses phase-referenced, and there- 
fore position-referenced, model components. If the target was 
not visible after phase referencing, the data were self-calibrated 
using a point source as the initial model. In this case, absolute 
position information was lost (see Table 2). 

Absolute source positions were measured from the phase- 
referenced, self-calibrated images, where available; other- 
wise, relative source positions were measured from the self- 
calibrated images. Positions were always measured by fitting a 
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two-dimensional parabola to the source and recording the po- 
sition of the peak. This was found superior to fitting Gaussians 
because the images of the sources were not always well de- 
scribed by a Gaussian. Phase-referenced positions are based 
on the calibrator positions in Table 1. Absolute position un- 
certainties are typically 10"^ times the source-calibrator sep- 
aration (Lestrade 1991). Our random position uncertainties 
should be a small fraction of the beamwidth, and proportional 
to beamwidth/SNR, where SNR is the signal to noise ratio of 
the source component. 

Flux densities were measured by integrating over the source 
region in the image plane. The estimated 1 cr flux density errors 
comprise a 5 % (MERLIN) or 10 % (other arrays) flux density 
scale calibration uncertainty, » 0.2 mJy additive thermal noise, 
the exact value depending on the source size and the particular 
image, and 5 % to 10 % uncertainty in the integration over 
the source region, depending on the size and flux density. The 
combined estimated errors of the integrated flux densities in 
Table 2 represent image rms noise plus 1 1 % (5 y < 20 mJy) and 
7 % (^v > 20 mJy) for MERLIN observations and 15 % (^v < 
20 mJy) and 1 1 % (5 y > 20 mJy) for other VLBI observations. 
The peak flux density errors are image rms noise plus 5 % for 
MERLIN observations and 10 % for other VLBI observations. 

Spectral indices were measured by tapering the 6 cm data 
and restoring with the 1 8 cm beam size. The uncertainty on the 
spectral indices, based on the uncertainties on the flux density 
measurements is 0.09 (1 cr) for a 20 mJy component. Despite 
the taper, some beam mismatch remained because the 18 cm 
and tapered 6 cm arrays were not perfectly scaled arrays; the 
resulting magnitude of the efifect is difficult to quantify because 
it depends on the source structure. The largest eff"ect probably 
comes from the short baselines at 18 cm, since these were not 
removed to match the 6 cm short-baseline length. This would 
aff"ect sources that have emission that is extended on scales 
larger than the largest angular scale to which the 6 cm ob- 
servations are sensitive (11 mas for the EVN), but still small 
enough to be detected at 18 cm (< 35 mas for the EVN). The 
insensitivity of the 6 cm observations to emission on this range 
of scales makes the spectral index appear more negative than 
it actually is. We see emission extended on this scale in the 
case of NGC 5506 (components BO and Bl) and NGC 7674 
(Momjian et al. 2003), and we attempt in these cases to esti- 
mate a correction. For the other sources our EVN 18 cm im- 
ages and MERLIN images do not show structure extended on 
this scale and so we beUeve that no bias in the spectral indices 
should be present. 

3. Results and Discussion 

3.1. NGC 5506 

NGC 5506 is an edge-on, dusty irregular or early-type spiral. 
It has formerly been classified as a type 2 Seyfert based on the 
nuclear optical spectrum (Wilson et al. 1976), but recent near- 
IR spectra suggest that it is a narrow-Une Seyfert 1 (Nagar et al. 
2002b). 

Colbert et al. (1996) detected large-scale radio emission ex- 
tending up to 3.6 kpc away from the disk in a north-south di- 



rection (i.e. perpendicular to the disk), but could not establish 
whether this is starburst or AGN driven. At higher resolution 
(~ 0.7"), a MERLIN 18cmimagebyUngeretal. (1986) shows 
that NGC 5506 has a compact core with a diffuse halo of a few 
arcseconds diameter. The compact core is slightly extended and 
aligns with the major axis of the host galaxy and with the posi- 
tion angle of the optical continuum polarization (Martin et al. 
1983). This is contrary to the situation in most Seyfert 2s, in 
which the polarization tends to be perpendicular to the jet axis 
(Antonucci & Miller 1985) and strengthens the Sy 1 classifica- 
tion. 

A water maser was discovered in NGC 5506 by Braatz 
et al. (1994) during a survey of Seyfert galaxies. Unfortunately, 
its weakness, 60 mJy in the broad component, has prevented 
detailed VLBI imaging. 



3.1.1. Results 




50 40 30 20 10 -10 -20 

miltiarcsec 



Fig. 1. VLBI 6 cm images at four epochs (left panels), from which we 
placed an upper limit on the proper motion of the Bl component, rel- 
ative to BO, in NGC 5506. The proper motion of the B2 component 
is unreliable due to limited image fidelity. The 6 cm 1994.16 epoch is 
repeated as the top right panel. The lower three right panels show the 
18 cm images. All images are naturally weighted. The VLBA obser- 
vations were restored with common beams of 7.4 x 4.3 mas (6 cm) and 
18.5x11.3 mas (18 cm). 

Both the MERLIN 18 cm and 6 cm images (not shown) dis- 
play a single unresolved point source with a total flux density 
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Fig. 2. Diagram illustrating the apparent motions in NGC 5506. The 
peak positions of B 1 and B2 from 6 cm global VLBI, EVN and VLBA 
are plotted relative to BO. Errors are 1 cr. 



of 1 13 mJy and 86 mJy, respectively. The single-dish 6 cm flux 
density in a 2.4' beam measured simultaneously at Effelsberg 
was 193 mJy, so fliere is significant emission on scales larger 
than the MERLIN 6 cm beam size, almost all of which appears 
in the VLA image of Colbert et al. (1996). 

The VLBI 18 cm and 6 cm images in Fig. 1 show fliree 
components which we name BO through B2, with separations 
of 3.5 pc (BO-Bl) and 5. 1 pc (B0-B2). These three components 
are unresolved or sHghdy resolved by the VLBA and EVN. The 
total 18 cm flux density recovered by the EVN in 1999.86 was 
65 mJy, one third of the 6 cm single-dish emission. 

We have aligned images from all four epochs using the cen- 
ter of BO as a reference point to illustrate the relative motions 
between components inside NGC 5506 (Fig. 2). From 1994. 16 
to 1997.10, component Bl has moved (1.60 ± 1.26) mas (1 cr 
error), yielding (0.21 +0.17) c. However, the error is large, and 
motion as fast as 0.50 c would be required for a 3 cr detection. 
After 1997.10, the motion of Bl is uncertain. The position an- 
gle in which Bl might be moving is (18+2)°, different from the 
P.A. of (72 ± 2)° of B 1 with respect to BO, so if the motion of 
Bl is real, it is not simply radially away from BO. Component 
B2 showed no consistent believable displacements; the appar- 
ent motions of B2 are probably due to Umited image fidelity at 
low SNR. 

The diameters of components BO and Bl were measured 
by fitting a Gaussian to the images, and the intrinsic size 
was determined by deconvolving the beam from the fitted 
Gaussian. Component BO proved to be only slightly extended 
(3.2 X 2.6 mas) at 6 cm compared to the beam size of 5.7 x 
3.7 mas (VLBA, February 1997), consistent with the near 
equality of peak and integrated flux densities. Component B 1 
was significantly resolved in the same observation (12.5 x 
8.9 mas), confirmed by the significantly smaller peak flux den- 
sity (4.3 mJybeam"') compared to the integrated flux density 
(24.4 mJy). The brightness temperature of component BO, de- 
rived from the integrated flux density and deconvolved compo- 



nent size, is 3.6 x 10** K, whereas the brightness temperature of 
the resolved component B 1 , derived from the peak flux density 
and beam size, is 1.4 x 10^ K. 

Spectral indices were measured using EVN 18 cm and EVN 
6 cm images tapered to matching resolution. The strongest 
component (BO) was found to have a flat spectrum (or = -0.03) 
while the other components had steeper spectra (a = -0.77 and 
-0.3). 

The spectral index measured for component BO might be 
affected by the extended emission to its NW (Fig 3, adding a 
negative bias to our measurement. To estimate the size of the 
effect, we note that the extended emission has a brightness of 
2mJybeam"' at 6 cm in a 16.5x8.0 mas beam, and is present in 
one quadrant from BO. This emission would contribute a flux 
density of 3 mJy in the 26.2 x 14.1 mas tapered beam area 
at 18 cm assuming it has a spectral index of -0.7 and that it 
would be seen by the EVN 18 cm array but not by the 6 cm 
array. Subtracting 3 mJy from the measured flux density of BO 
at 18 cm of 32.9 mJy yields a corrected spectral index of +0.06, 
as listed in Table 2. 

Likewise, the spectral index estimated for component Bl 
might be affected by a possible bridge of faint emission be- 
tween Bl and BO visible in the combined EVN and MERLIN 
image in Fig. 3. The emission has a brightness of 2 mJy beam"' 
at 6 cm and occupies one octant. This emission would con- 
tribute 1.5 mJy to the EVN 18 cm measurement and would be 
resolved out by the 6 cm EVN array. Subtracting these 1 .5 mJy 
from the measured flux density of Bl at 18 cm of 25.6 mJy 
yields a corrected spectral index of -0.73 (Table 2). 

The flux density of component BO remained constant to 
within the uncertainties. Component B 1 faded by a factor of 
two at 6 cm and changed its spectral index from a^^ = 0.0 in 
1997 to = -0.77 in 2000. The change could be due to an ex- 
panding synchrotron source that was initially compact enough 
to be synchrotron self-absorbed, or else could be due to chang- 
ing foreground free-free absorption. 

Combining the EVN and MERLIN 6 cm data yields a pos- 
sible bridge of emission connecting the two brightest compo- 
nents and the detection of very faint diffuse structure north of 
component BO (Fig. 3). This diffuse structure is interesting be- 
cause it suggests a physical connection between the AGN and 
the radio emission on kpc scales observed with the VLA in 
C-configuration by Colbert et al. (1996). The MERLIN 18 cm 
beam size of 260 x 150 mas includes the three components BO- 
B2 and the diffuse northern emission (80 X 40 mas) within a 
single beam. Thus, the total MERLIN 18 cm flux of 1 13 mJy 
is probably the combined emission from all the structures seen 
at 6 cm. This diffuse northern emission is not seen in the EVN 
image because the EVN is insensitive to structures larger than 
35 mas at 18 cm or 11 mas at 6 cm. Although component 
B2 has the same brightness as the northern diffuse region in 
the MERLIN image, B2 remains visible in the EVN images 
whereas the northern diffuse emission does not, because B2 is 
more compact and so is not resolved away at the higher EVN 
resolution. 
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Fig. 3. Combined EVN and MERLIN 6 cm image of NGC 5506. 
Offsets are relative to RA 14:13:14.87953, DEC -03:12:27.6799 
(J2000), peak flux density is 37.3 mjy/beam, contours are 
1 mjy/beamx(2, 4, 8, 16, 32), beam size is 16.5x8.0 mas and natural 
weighting was used. 



3.1 .2. Discussion 

The nuclear radio emission of NGC 5506 is resolved into three 
compact, non-colinear core components within the central 
5 pc, and a diffuse region of 8 pc extent to the north. The 
two strongest radio components (BO and Bl) are misaligned 
by 20° from the galactic plane and almost perpendicular to 
the kpc radio structure. On the hundred pc to kpc scale, the 
radio jet axes in Seyferts have been found to be randomly 
oriented with respect to the galactic disk rotation axes (e.g. 
Pringle et al. 1999, Nagar & Wilson 1999, Kinney et al. 
2000), whereas by conservation of angular momentum during 
accretion, one would expect a common axis. Misaligrmient of 
pc and kpc-scale radio emission is shown later to be a common 
phenomenon in Seyfert galaxies (Table 8, Fig. 7(b)). 

The radio emission from component B lis peculiar. Its 6 cm 
integrated flux density decreased from 30.0 mJy in 1994 to 
24.4 mJy in 1997 and 10.8 mJy in February 2000, and flien 
increased slightly to 15.3 mly in May 2000 (this last increase, 
however, is not significant within 3 a). Thus, the 6 cm flux den- 
sity has constantly decreased, producing a corresponding steep- 
ening of the spectrum over this six year period from a = 0.0 
in February 1997 to a = -0.78 in May 2000. The combina- 
tion of flat spectrum and low brightness temperature in 1997 
(1.4 X 10^ K) cannot be explained by synchrotron self absorp- 
tion, and is probably due to free-free absorption. The time vari- 
abiUty of the spectrum rules out free-free emission because the 



spectral cutoff of a gas with such a high temperature is in the 
X-ray regime and radio spectral index variability as the com- 
ponent cools is therefore not expected. 

BO is a candidate for the location of the AGN core due to 
its flat spectrum (a^^ = 0.06) and compactness. 

The alternative picture is that the AGN might lie between 
BO and B 1 and then the two bright components might be work- 
ing surfaces where double-sided jets from the AGN impact the 
ISM. Low-surface -brightness emission is seen extending from 
the hot spots to larger scales (the region of emission north of 
BO, and component B2) and could be plasma that has flowed 
through the hot spots and is either deflected by the shock, or is 
rising buoyantly in the pressure gradient in the ISM. However, 
BO has a flat radio spectrum, indicative of optically thick in- 
stead of optically thin synchrotron emission, whereas working 
surfaces, like in compact symmetric objects (CSOs), charac- 
teristically show steep spectra. Hence, we favour BO being the 
core. 

The picture of an AGN core lying between two hotspots 
seems to describe the radio observations of NGC 7674 dis- 
cussed in the next section and in Momjian et al. (2003), and 
of NGC 3079 (Middelberg et al., in prep.). 

BaUistic motion of components, when back-extrapolated, 
should converge on the core. However, in NGC 5506 the Bl 
component, if its motion is real, moves almost northwards in 
P.A. (18 + 2)°, whereas its position with respect to the likely 
core component BO is in P.A. (72 + 2)°. 

A slight variant on this picture would have precession of the 
accretion disk move the hot spots "sideways" along a "wall" 
of ISM leaving behind a trail of fading extended plasma to ac- 
count for the observed proper motions having a tangential com- 
ponent, and for the "trail" of low surface-brightness emission 
behind each hot spot. 

3.2. NGC 7674 

NGC 7674 is a face-on spiral galaxy with asymmetric arms 
and a tidal connection to a nearby compact elliptical galaxy 
(Hutchings 1996). Observations witii the VLA and EVN 
(Unger et al. 1988) reveal a linear triple radio source of ~ 0.7" 
angular extent. The components are ~ 0.5" west and ~ 0.15" 
east of the main peak. Momjian et al. (2003) made a very sen- 
sitive observation of NGC 7674 with the VLBA, phased VLA 
and Arecibo at 1.4 GHz showing the triple source found by 
Unger et al. (1988) and additional low surface-brightness emis- 
sion forming an S-shaped structure. 

3.2.1. Results 

Our MERLIN 18 cm image (Fig. 4(a)) shows the well-known 
triple structure consistent with the VLA A-array 15 GHz im- 
age by Unger et al. (1988), but the eastern component was not 
detected by the EVN or by MERLIN at 6 cm. 

The integrated MERLIN 6 cm flux density was 25.6 mJy, 
lower than the total 6 cm single-dish flux density of 73 mJy 
in a 2.4' beam. The latter was measured three months earlier 
at EfFelsberg in a related project to look for Seyfert variabil- 
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(a) NGC 7674 at 18 cm (MERLIN). Offsets are relative to 
RA 23:27:56.71696, DEC 08:46:44.0078 (J2000), peak flux den- 
sity is 74.9 mjy/beam, contours are I mJy/beamx(-5, 5, 10, 15, 
20, 75 ), beam size is 276.1 masx 110.4 mas and uniform 
weighting was used. 



milliarcsec 
-200 -300 




-200 -300 
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(b) EVN 18 cm images of NGC 7674. Top: First epoch at 1985.28. Peak flux den- 
sity is 29.8 mJy/beam, contotxrs are 1 mJy/beamx(-2, 2, 4, 8, 16, 25), beam size 
is 39.3 masx24.2 mas and natural weighting was used. Bottom: Second epoch at 
1999.86. Offsets are relative to RA 23:27:56.71201, DEC 08:46:44:1368 (J2000), 
peak flux density is 23.9 mJy/beam, contours are 1 mJy/beamx{— 1, 1, 2, 4, 8, 16, 
32), beam size is 32.5 masx9.8 mas and uniform weighting was used. 



Fig. 4. NGC 7674 images 



ity. Thus, only 35 % of the total emission was detected by the 
MERLIN observation; the remaining 65 % is also not seen in 
the VLA 15 GHz image of Unger et al. (1988). Both MERLIN 
and VLA observations are insensitive to structures larger than 
0.15 arcsec, or 90 pc. 

Our 18 cm EVN image from November 11, 1999, in 
Fig. 4(b) (lower panel) displays a bright, slightly resolved com- 
ponent and point-like radio emission towards the north west, 
consistent with the structure seen previously by Unger et al. 
(1988). The SE component is slightly extended in the direction 
of the NW component. 

From the 1985 EVN archive observation we obtained a 
good-quality 18 cm image (Fig. 4(b), upper panel) with a 49 cr 



detection of the SE component and a 9 cr detection of the NW 
component. We produced images from the 1985 and the 1999 
data sets with matching (m, v) coverages and equal restoring 
beams (41 .9x24.5 mas^ in RA. -29.20°) to look for a change in 
separation of the components. From the 1985 image, we mea- 
sured a separation of (508.7 +0.4) mas (1 cr error), and from the 
1999 image, we measured (516.0 + 0.4) mas (1 cr error). Thus, 
we find that the relative separation between the components has 
increased by (7.3 + 0.6) mas in 14.58 years, corresponding to 
0.28 pc yr ' or Vapp = (0.92 ± 0.07) c. The separation increased 
along P. A. 110°, in agreement with the P. A. of the the cen- 
tral and the NW component of 1 17°. The 18 cm flux densities 
of the SE and NW component have increased from 22.3 mJy 
to 37.6 mJy and from 4.1 mJy to 11.0 mJy, respectively, the 
second epoch being in agreement with Momjian et al. (2003). 
Amplitude calibration errors could in principle produce a scal- 
ing of the apparent flux densities, but probably not by such a 
large factor as seen here. Furthermore, they cannot produce the 
observed decrease of the ratio of the flux densities from 5.4 to 
3.4. Thus, it is most likely that both components have bright- 
ened. 

The spectral indices were calculated using the EVN 18 cm 
and EVN 6 cm images tapered to matching resolution. The low 
surface brightness emission surrounding component NW seen 
by Momjian et al. (2003) (1.7 //Jy mas^) has a neghgibly small 
effect on the spectral indices measured by us. 

3.2.2. Discussion 

All components in the EVN and MERLIN images of 
NGC 7674 have steep spectra, and no core identification can 
be made based on the spectrum. Momjian et al. (2003) suggest 
that the nucleus might lie between the two components we de- 
tected and is free-free absorbed at 2 1 cm. If so, our observations 
show that it must also be absorbed at 6 cm. 

If the core is indeed located between the two brightest 
components, then the structure would be reminiscent of that 
in compact symmetric objects, and the components could be 
working surfaces where the jets impact the ISM. Then the 
proper motions would measure the hot-spot advance speed. 
From the proper motion of 0.14 pcyr"^ and the component 
separation of 282 pc we derive a dynamical age for the source 
of 2000 yr, which lies within the range of ages measured for 
CSOs by Polatidis et al. (2002). However, the 18 cm radio 
luminosity of the two components in NGC 7674 combined 
is 8 X 10^^ Wm^^Hz"', which is three orders of magnitude 
less than the typical CSOs studied by Polatidis et al. (2002). 
The high radio luminosity of CSOs is probably a selection 
effect, since the CSO samples are drawn from flux-density- 
limited surveys and so automatically exclude lower-luminosity 
sources. Thus, we may here be observing the low-luminosity 
end of the CSO phenomenon. 

The speed of (0.92 ± 0.08) c with which the two compo- 
nents separate is, after inZw2, the highest found in a Seyfert 
galaxy. If the components are shocks across jets, then the flow 
speed of the jet could be even higher. Furthermore, the PA. in 
which the separation increases is in agreement with the relative 
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P.A. of the central and the NW component, indicating a linear 
expansion of the source. The separation of 282 pc between the 
components is comparable to the longest jets found in Seyferts 
(e.g. NGC 4151, Pedlar et al. 1993 or Mrk 3, Kukula et al. 
1993), and shows that jets in Seyfert galaxies can propagate 
large distances while maintaining high speed. 

NGC 7674 is remarkable for its broad asymmetric for- 
bidden line profiles which show highly blue-shifted wings 
(2000 kms"', Unger et al. 1988). These blue wings lie well 
outside the galaxy rotation curve and a connection to the AGN 
or jet seems Ukely (Unger et al. 1988), as has been found 
in other objects (e.g. Whittle et al. 1988, Axon et al. 1998). 
Whittle (1992) investigated the acceleration mechanisms of 
NLR gas in Seyfert galaxies. In his sample, Seyferts such as 
NGC 7674 with linear radio morphology and radio powers ex- 
ceeding 1022-5 w Hz-i (NGC 7674 has 2.5 x lO^^ W Hz ', de- 
rived from MERLIN 18 cm integrated flux density) have sys- 
tematically broader [O III] lines, which he attributes to jet ma- 
terial acting on the NLR gas. 

3.3. NGC 2110 

NGC 2110 is an elliptical or SO galaxy. VLA observations at 
6 cm and 20 cm with 0.4" and 1.2" resolution by Ulvestad 
& Wilson (1983) show a symmetric, S-shaped radio source 
~ 4" (600 pc) in extent which was interpreted as a two-sided, 
continuous jet interacting with the ISM. This structure was 
confirmed in a VLA 3.6 cm image by Nagar & Wilson (1999) 
which also resolved finer details in the jets. 



3.3.1. Results 

Our MERLIN 18 cm image (Fig. 5, tapered to 0.3 x 0.2 arcsec 
resolution shows the core and two components, to north and 
south, at distances of (1.72 ± 0.02)" and (1.62 ± 0.02)", re- 
spectively, from the central component. These separations are 
consistent with those measured by Ulvestad & Wilson (1983) 
for the brightest parts of the radio lobes. The extended struc- 
ture seen in their observations was not detected by EVN 
or MERLIN as the surface brightness (360 K) is below the 
MERLIN detection threshold (2000 K), except for the bright 
hot spots. The MERLIN 6 cm observations yielded only detec- 
tion of a point source (not shown). 

The total flux density in the MERLIN 6 cm image was 
37 mJy, compared to the single-dish 6 cm flux density of 
154 mJy in a 2.4' beam measured simultaneously at Effelsberg. 
Thus, three quarters of the total flux density originates from 
emission extended over more than 0.7" or 75 pc. The total VLA 
6 cm flux density was (175 ± 8) mJy in 1983 with a 0.4" beam 
(Ulvestad & Wilson 1983), which is marginally more than the 
single-dish flux density in 2000. 

The full-resolution NGC 21 10 EVN images at 18 cm and 
6 cm (not shown) display only an unresolved point source at 
the position of the VLA observed core. The spectral index of 
the core was measured using an EVN 6 cm image tapered to 
a resolution matching that of the EVN 18 cm image and the 




arcsec 



Fig. 5. NGC 2110 at 18 cm (MERLIN). Offsets are relative to 
RA 05:52:11.38254, DEC -07:27:22.5262 (J2000), peak flux density 
is 31.4 mjy/beam, contours are 1 mJy/beamx(-3, 3, 6, 9, 30), 
beam size is 318.0 masxl93.4 mas and natural weighting was used. 
The (m, v) range was 80 - 1000 kA. 

spectrum was found to be rising slightly with increasing fre- 
quency (a - +0.12). The MERLIN spectral index of the core 
(-hO.20) was calculated using the full (u, v) range because ta- 
pered 6 cm images had too little {u, v) coverage remaining to 
yield reliable images, and no detection was made of the north 
and south components. 

Since there is no extended emission on the 11 mas to 
35 mas scale near the core, we believe that these spectral in- 
dices are unbiased by the short-basehne-mismatch between 
the arrays. Observations at higher resolution carried out by 
Mundell et al. (2000) with the VLBA at 3.6 cm confirm that 
the core is unresolved down to a resolution of 2 mas (0.3 pc). 

3.3.2. Discussion 

The spectral index of the core in NGC 2110 was estimated by 
Ulvestad & Wilson (1983) and Nagar et al. (1999) from VLA 
observations, giving = 0.0 and aj g = -0.47, respectively. 
However, both measurements suffer from the low resolution of 
the VLA at 20 cm, which does not allow a clean separation of 
core and lobes. Therefore, our EVN and MERLIN observations 
have aUowed an accurate measurement of the core spectral in- 
dex for the first time. 

3.4. Mrk 1210 

Mrk 1210 is an almost face-on, amorphous Seyfert 2 galaxy. 
It has a compact NLR (Falcke et al. 1998), and shows broad 
emission lines in polarized light, indicating a hidden type 1 
AGN (Tran et al. 1992). Mrk 1210 has also been classified as 
a Wolf-Rayet-Galaxy due to its bright Hell /t4686 emission 
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(Storchi-Bergmann et al. 1998). Falcke et al. (1998) imaged 
Mrk 1210 in Hot and [Om] with the HST, to try to resolve an 
ionization cone, and at 3.6 cm with the VLA, with resolutions 
of 100 mas and 200 mas, respectively. In both of those obser- 
vations Mrk 1210 appears compact though with a marginally 
significant extension towards the southeast in both their optical 
and radio observations. 

The single-dish 6 cm flux density was 76 mJy in a 4.2' 
beam measured in 1990 (Griflith et al. 1995) at the Parkes 
64 m telescope. Mrk 1210 hosts an H2O megamaser with an 
isotropic luminosity of 84 Lq (Braatz et al. 1994) which might 
indicate the presence of an edge-on circumnuclear molecular 
disk or torus. 



3.4.1. Results 

The MERLIN 18 cm image (not shown) displays an unresolved 
point source with a possible weak halo to the north-west. The 
MERLIN 6 cm image (Fig. 6(a)) reveals a bright, unresolved 
central component with a weak extension to the north-west, 
and a third component to the far south-east, which we name F. 
The total MERLIN 6 cm flux density was 39 mJy compared 
to tiie 6 cm single-dish flux density of 54 mJy in a 2.4' beam 
measured simultaneously at Eff'elsberg on the same day, and so 
the MERLIN image has recovered most of the flux density. 

The EVN and VLBA 18 cm images in Fig. 6(b) and 6(c) 
reveal a bright compact object (named NW) and a weaker 
component towards the south-east (named SE) at a distance of 
(32.7 + 1.0) mas (8.6 pc) from the NW object. We marginally 
detect a third component at a distance of » 117 mas (30.6 pc) 
from NW which is spatially coincident with the weak object 
F in the MERLIN 6 cm image. This component is not for- 
mally significant but is present in both images and so is proba- 
bly real and could be the continuation of the radio ejecta. The 
EVN 6 cm data (not shown) yield only a marginal detection of 
the strongest point source. The 6 cm VLBA image (Fig. 6(d)) 
resolves the brightest source into an arc of four components 
and also resolves the bright component, named SE. The total 
flux density recovered by the full-resolution EVN and VLBA 
at 18 cm was 50 mJy and 49 mJy, and at 6 cm was 10 mJy and 
9.3 mJy, respectively. The total 6 cm flux density measured by 
the VLBA in a beam tapered to the 18 cm resolution is 14 mJy, 
which is stiU only 35 % of the total MERLIN flux density. So 
most of the galaxy's 6 cm emission comes from scales between 
3 pc (the largest angular scale to which the EVN is sensitive) 
and merlin's beam size of 16 pc, and only little from even 
larger scales. 

Both the NW and SE components have steep spectra, with 
a = -1.26 and a = -0.78 in our VLBA images. In our EVN 
images, the NW component has a = -1.19 provided that the 
weak detection in our EVN 6 cm observations is the NW com- 
ponent. Since aU components have steep spectral indices, the 
core identiflcation is best done based on proper motions. Our 
observations here provide the first epoch absolute positions and 
a second epoch should follow after at least five years to allow a 
reliable detection of proper motion. 



The position angle of the SE component with respect to 
the NW component (122 ± 5°) is perpendicular to the optical 
polarization vector (RA. 29°) measured by Tran (1995), as is 
generally the case in Seyfert 2s. 

3.4.2. Discussion 

We suggest three scenarios for the radio structure in Mrk 1210. 

A free-free emitting disc? The component configuration 
reminds us of NGC 1068 in which the core component SI at 
the jet base is resolved into a 15 mas (1 pc)-long structure per- 
pendicular to the jet axis in the sensitive 3.6 cm VLBA obser- 
vation by Gallimore et al. (1997). Surrounding this component 
is an edge-on disc of H2O masers whose rotation curve implies 
a compact central mass of 1.5 x 10^ Mq (Greenhill & Gwinn 
1997). In that case, the 3.6 cm continuum emission was iden- 
tified as free-free emission from the ionized inner edge of the 
torus, based on its flat spectrum, relatively large extent, orien- 
tation perpendicular to the jet, and low brightness temperature 
(10^ K), the last ruling out synchrotron self-absorption as the 
cause of the flat spectrum. In Mrk 1210, the bright arc of four 
components suggests an edge-on disk, whose existence is fur- 
ther supported by the presence of nuclear water maser emission 
(although the precise location of the masers is not known). The 
integrated radio spectrum of component NW is a^^^ - -1.2 us- 
ing the EVN and a^^ = -1.3 using the VLBA. This favours 
optically thin synchrotron radiation rather than flat-spectrum 
free-free emission. The brightness temperature measured with 
the VLBA at 18 cm is 9.2 x 10^ K and with the EVN at 18 cm 
is 7.1 X lO'' K, significantly higher than in NGC 1068. Thus, 
although the shape of the radio emitting regions in component 
NW is suggestive, the structure cannot be free-free emission 
from the inner edge of a torus like in NGC 1068. 

Synchrotron emission from a torus or extended accre- 
tion disk? The arc of four components in the NW compo- 
nent could be synchrotron emission from an edge-on viewed 
torus or accretion disc. However, accretion disks are normally 
thought to be much smaller and neither accretion disks nor tori 
are known as synchrotron emitters. 

A core and outflow? Perhaps component SE hosts the nu- 
cleus. The resolved arc of four components in NW then could 
be bright knots due to a shock in an extended, low brightness 
temperature radio outflow, with similar conditions to those in 
NGC 5506, NGC 1068 or Mrk 231. But both the SE and flie F 
component lack a flat spectrum that would indicate the nucleus. 

Further clues might come from absolute proper motion 
measurements in the future, or from VLBI imaging of the water 
maser emission. 

4. Statistics 

The range of radio properties displayed by Seyfert galaxies has 
been explored extensively based on WSRT, VLA or ATCA ob- 
servations at 20 cm, 6 cm or 3.6 cm at arcsec resolution by 
de Bruyn & Wilson (1978), Ulvestad & Wilson (1984, 1989) 
Unger et al. (1987), Giuricin et al. (1990), Rush et al. (1996), 
Nagar & Wilson (1999), Morganti et al. (1999), Ulvestad & Ho 
(2001a), and Vir Lai (2001). They studied the range of radio 
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(a) MERLIN 6 cm image of Mrk 1210. Offsets are relative to 
RA 08:04:05.8601 1, DEC 05:06:49.8992 (J2000), peak flux density 
is 24.1 mjy/beam, contours are 1 mjy/beamx(-l, 1, 3, 5, 7, . . ., 23), 
beam size is 72.7 masx47.6 mas and uniform weighting was used. 



(b) EVN 18 cm image of Mrk 1210. Ofl'sets are relative to 
RA 08:04:05.856, DEC 05:06:49.83 (J2000), peak flux density is 
19.8 mJy/beam, contours are 1 mjy/beamx(-l, 1, 3, 5, 19), 
beam size is 24.6 masxl 1.3 mas and uniform weighting was used. 
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(c) VLBA 18 cm image of Mrk 1210. Offsets are relative to 
RA 08:04:05.841, DEC 05:06:49.72 (J2000), peak flux density is 
21.6 mJy/beam, contours are 1 mJy/beamx(-l, 1, 3, 5, .... 21), 
beam size is 18.6 masxl4.2 mas and natural weighting was used. 



(d) VLBA 6 cm image of Mrk 1210. Offsets are relative to 
RA 08:04:05.841, DEC 05:06:49.72 (J2000), peak flux density is 
2.5 mJy/beam, contours are 0. 1 mJy/beamx(-3, 3,6,..., 24), beam 
size is 3.2 masxl. 5 mas and uniform weighting was used. 



Fig. 6. Mrk 1210 images. 



luminosities, radio source sizes, morphologies, alignment be- 
tween jet and host galaxy, and spectral indices, and looked for 
dependence of properties on Seyfert type. With the exception 
of Ulvestad & Ho (2001a), who found 50% of low-luminosity 
Seyferts to have a flat spectrum core, they all found that most 
Seyfert galaxies have steep, optically-thin synchrotron radio 
spectra, with 5 % to 35 % (depending on the sample) showing 
a flat-spectrum core, much less than the rate of occurrence of 
flat-spectrum cores in powerful radio galaxies. No difference 
was found between the spectral index distributions of the two 
Seyfert types. 



However, the 0.3" to several arcsecs beams used in those 
studies do not allow details of nuclear structure to be sepa- 
rated or the distribution of jet bending angles to be constructed. 
Further, they would probably blend together multiple radio 
components within the nucleus and could hide the presence of 
more absorbed cores. Our VLBI observations provide enough 
resolution to measure separately the emission from individ- 
ual sub-pc-scale components, and so can find parsec-scale ab- 
sorbers more efiiciently. However, our sample size is too smaU 
for any general conclusions. 

Many further VLBI observations of Seyfert galaxies are 
available in the literature but no statistical summary from those 
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observations has been presented. Thus, we made a literature 
search for VLBI observations of Seyferts to supplement our 
sample. The resulting list in Table 4 contains, we believe, all 
VLBI imaging observations of AGN classified as Seyferts by 
Veron-Cetty & Veron (2001), pubUshed as of August 2003 
with a linear resolution of < 10 pc, although it is not impossi- 
ble that we have inadvertently omitted some works. This sam- 
ple consists mostly of the nearest and best-known Seyferts. 
They were often selected for observing because they show 
bright radio nuclei at arcsecond resolution, and so the sample is 
not complete in any sense. Large, uniform samples of Seyfert 
galaxies are being observed with the VLBA by Murray et al. 
(1999) and Schmitt et al. (2002), but the results are not yet 
available. Proper motions in Seyferts are being presented by 
Roy et al. (in prep.) and we confine ourselves here to spectral 
indices and jet bending. 



2 

Spectral Index 



4.1. Fraction of Flat-Spectrum Components 

We selected from Table 4 a subsample consisting of the Seyfert 
galaxies for which dual-frequency VLBI observations have 
been made, including the four objects presented in this paper, 
to study the component spectral properties at high resolution. 
This subsample is listed in Table 5. 

Of the 16 objects in the subsample, a majority of 13 (81 %) 
has at least one flat or inverted spectrum VLBI component with 
ttg^ > -0.3. This fraction is 60% larger than that measured 
at lower resolution by Ho & Ulvestad (2001) for the Palomar 
Seyferts and four or more times larger than the fraction mea- 
sured at arc sec resolution for the other samples (Table 6). The 
diff^erence is highly significant for most samples, having a prob- 
ability of < 0.1 % of occurring at random. The larger fraction 
of flat- and inverted-spectrum components found with VLBI is 
not unexpected since increasing resolution reduces the compo- 
nent blending. 

4.2. Spectral Index Distributions 

The distribution of a^^ from Table 5 is shown in Fig. 7(a). 
The spectral indices lie between -1.8 and +4.6 with a me- 
dian of -0.40 (using the Kaplan-Meier estimator constructed 
by ASURV rev 1.2, LavaUey et al. 1992 to account for lower 
limits and discarding upper limits). The distribution shows a 
tail towards positive spectral index due to four components 
in NGC 1068, NGC 2639 and NGC 3079 tiiat have strong 
low-frequency absorption. In the case of two components in 
NGC 3079, the spectral indices are greater than the maximum 
a^^ = +2.5 that can be produced by synchrotron self absorp- 
tion and so is most likely due to free-free absorption, probably 
by dense circumnuclear plasma ionized by the AGN. 

The tail of extreme positive spectral indices that we see in 
our sample is not seen at arcsec resolution in the works men- 
tioned above, probably due to component blending. We made 
a statistical comparison of our VLBI spectral indices to arc- 
sec resolution spectral indices using the Kolmogorov-Smimov 
test and chi-squared homogeneity test by treating the lower and 
upper limits as detections, and by applying survival analysis 
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Fig. 7. Upper panel: Histogram of spectral indices between 18 cm and 
6 cm for the Seyfert sample in Table 5. Upper and lower limits are 
indicated by arrows. Lower panel: Distribution of angles between pc- 
and kpc-scale structure. For comparison with data from radio galaxies, 
the data from Table 7 have been reflected onto the range 0° to 90°. 



two- sample tests. The results are summarized in Table 7. Our 
spectral index distribution was significantly difi'erent (i.e. level 
of significance was less than the critical value of 5 %) from 
all the comparison samples according to all tests performed, 
with the exception of the Palomar Seyferts by Ho & Ulvestad 
(2001), which was found to be significantly diff'erent in only 
the chi-squared test. This excess of flat-spectrum sources in the 
Palomar Seyfert sample was already noticed by Ho & Ulvestad 
(2001). 
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Separating the VLBI measurements by Seyfert type, we 
found the median a^^ for the Sy Is was -0.75 and for the Sy 
2s was -0.23, which were not significantly different under the 
survival analysis two-sample tests (f = 6 % to 12 % that the 
spectral-index distributions came from the same parent distri- 
bution) or under a Kolmogorov-Smirnov test (P = 6.1 % to 
come from the same parent distribution), treating only lower 
limits as detections. 

Despite the statistical insignificance of a difference between 
Seyfert types, we notice that the most inverted spectra all occur 
in objects that have evidence for an absorber. NGC 3079 is a 
type 2 Seyfert, NGC 1068 and Mrk 348 are type 2 Seyferts with 
type 1 spectra in polarized hght and NGC 2639 shows strong 
X-ray absorption (Wilson et al. 1998). The inner edge of such 
an absorber is expected to be ionized by UV emission from 
the AGN which would produce free-free absorption and hence 
produce inverted spectra preferentially in the type 2s. All four 
objects with inverted radio spectra also host H2O masers, which 
is also consistent with an orientation effect if, for example, the 
radio core is viewed through a molecular disk or torus, which 
produces H2O maser emission. 

4.3. Jet Misalignment 

Kiloparsec-scale linear radio structures in Seyfert galaxies are 
thought to be energized by jets or outflows originating at the 
AGN core, and so one might expect pc and kpc-scale radio 
structures to be aligned. Misalignments could be caused by 
changes in the ejection axis or by pressure gradients in the ISM. 
To look for such effects, we compiled kpc-scale radio observa- 
tions of the VLBI Seyfert sample, where available. The result- 
ing sample of 21 objects in Table 8 is, we beheve, a complete 
Ust of Seyferts with both pc-scale and kpc-scale observations 
published, although seven of these objects either have no clear 
extended structure on one of the scales or the structure was 
deemed to be affected by processes on scales > 1 kpc, e.g., 
buoyancy forces. In one object, NGC 1167, the double-sided 
kpc-scale emission is bent, so we quote both. This sample has 
a large overlap with the sample in Table 5 and, likewise, is not 
complete in any astronomical sense. Biases such as selection 
for bright radio nuclei are probably present, but whether they 
affect the distribution of jet bending angle is unclear. Doppler 
boosting would favour jets that point towards us, but the jet 
speeds so far measured in Seyferts are low and boosting effects 
should be minor. 

The distribution of misalignment angles from Table 8 is 
shown in Fig. 7(b). Misahgnment in Seyfert jets between pc 
and kpc-scales is common, with 5 out of 15 objects bending 
by 45° or more. The null hypothesis that the jet bending-angle 
distribution was drawn from a parent population of aligned jets 
was rejected by a chi-squared test at the 2 % level of signifi- 
cance. The null hypothesis that the distribution was drawn from 
a uniformly distributed parent population could not be rejected 
by a chi-squared test, which returned a 25 % level of signifi- 
cance. 

We compared the bending angles in Seyferts to those 
in radio-loud objects using a comparison sample of core- 



dominated radio-loud objects compiled from Table 6 in 
Pearson & Readhead (1988), and Table 1 in Conway & Murphy 
(1993), which have both VLBI and VLA or MERLIN observa- 
tions at 5 GHz. Those objects show a surprising bimodal dis- 
tribution of jet bending angle, in which the pc-scale jets tend 
to align with or to be perpendicular to the kpc-scale structure. 
We compared our jet bending angle distribution to the distri- 
bution in radio-loud objects using a Kohnogorov-Smimov test 
and found no significant difference, the probability of being 
drawn from the same parent distribution being 14 %. However, 
the Seyfert sample is small and the distribution for the Seyferts 
was shown in the previous paragraph to be consistent with a 
uniform distribution. 

The frequent misaUgnments seen in Seyfert galaxy jets 
could be evidence either for changes in the jet ejection an- 
gle and hence of the accretion disk, perhaps by a mechanism 
such as radiation-driven warping (Pringle 1997), or for bend- 
ing due to pressure gradients in the ISM (collision or buoyancy) 
through which the jets propagate. 

We compared the extended pc-scale radio structures to the 
host galaxy rotation axes, assuming that the host galaxy is cir- 
cular and the rotation axis vector is projected onto the galaxy's 
minor axis, which can be measured from optical images. For 
all spiral galaxies with pc-scale radio structure, we used the 
major axes from de Vaucouleurs et al. (1991), subtracted 90° 
and reflected the angles onto the range 0° to 90°. Four objects 
were not listed in de Vaucouleurs et al. (1991), so we used data 
from Schmitt & Kinney (2000) in two cases and we measured 
the major axes from Hubble Space Telescope archival data in 
the other two cases. We performed a Spearman rank correlation 
test and found a correlation coefficient of 0.07, corresponding 
to a 82 % significance and so there is no correlation between 
the directions of the pc-scale structures and the host galaxy ro- 
tation axes. However, the sample was small (13 objects) and 
therefore the test is not strong. 

5. Conclusions 

We have observed the Seyfert galaxies NGC 5506, NGC 7674, 
NGC 2110 and Mrk 1210 at 18 cm and 6 cm with the EVN 
and MERLIN. NGC 5506 and Mrk 1210 were additionally 
observed with the VLBA. These observations provided high 
angular resolution and good surface brightness sensitivity, 
spectral index measurements, and up to four epochs of obser- 
vations for proper motion measurements. Our conclusions are 
as follows. 

1. The structure seen on the tens-of-parsec-scale in 
NGC 7674 shows two bright, steep- spectrum components that 
might be the working surfaces where double-sided jets ejected 
from an (unseen) AGN core located between them impact 
the interstellar medium. Low-surface-brightness emission was 
seen extending from the hot spots in the direction of VLA- 
observed emission by Momjian et al. (2003). The diffuse emis- 
sion could be from material that has flowed through the hot- 
pots and was then either deflected or rose buoyantly in the 
pressure gradient in the ISM. Alternatively, the bend after the 
hot spot might be caused by accretion disk precession, which 
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might move the hot spots sideways along a wall of ISM leav- 
ing behind a trail of fading extended plasma. These structures 
are reminiscent of those seen in compact symmetric objects, so 
we suggest that NGC 7674 is a low-luminosity member of the 
CSO family. 

NGC 5506 shows two bright radio components separated 
by 3.5 pc with low-surface- brightness emission on larger scales 
extending towards the VLA-scale emission. The brightest com- 
pact component (BO) has a flat radio spectrum (or = 0.06) and 
could represent the nucleus. 

In NGC 21 10 and Mrk 1210, the overall structure is a triple 
radio source. For NGC 2110, the flat spectrum (a - 0.12) and 
the symmetric structure of the VLA-observed jets and the com- 
ponents at the ends of these jets (also seen in our MERLIN im- 
age) show conclusively that the central source represents the 
nucleus. In Mrk 1210, both the SE and the NW components 
have steep radio spectra and here the location of the nucleus is 
unclear. 

If jets in Seyferts are initially relativistic and are then 
disrupted by interaction with the ISM, we might expect 
significant faint, difiiise radio emission. Such faint emission 
(brighter than 2 mJy beam ') was seen in only NGC 5506. 

2. A proper motion of (0.92 + 0.07) c was measured in 
NGC 7674, and a 3 cr upper limit of 0.50 c was found for 
the motion of Bl in NGC 5506 with respect to BO. From 
the high speed measured in NGC 7674 between components 
282 pc apart we derived a dynamical age of 1000 yr. This 
speed and age is similar to those seen in CSOs. That the 
components have high speed after propagating hundreds of 
parsecs shows that the ISM in this object has not stopped 
the jet very effectively. The speed measured in NGC 7674 
is faster than those measured in other Seyferts by Ulvestad 
et al. (1998), Ulvestad et al. (1999b) and Roy et al. (2000), 
except for IIIZw2, though is slow compared to the speeds at 
equivalent scales in powerful radio galaxies. This suggests, but 
does not prove, that the jets are launched with an intrinsically 
lower speed than in radio galaxies. Further support comes from 
the preference of Seyfert nuclei for early-type spirals, which 
have relatively low-density ISM even close to the nucleus, and 
hence it seems unlikely that the jets are significantly slowed. 

3. We mentioned in the introduction the many possible 
reasons for low jet power in Seyferts that have been suggested 
in the literature. Our observations of proper motions favour 
an intrinsic over an extrinsic cause since the ISM has not 
stopped the component in NGC 7674, and no emission from 
disrupted jets was seen in three of four Seyferts. Any extrinsic 
deceleration must happen on scales less than 3.5 pc. (However, 
the bends in the radio sources at distances of 3.5 pc and 282 pc 
in NGC 5506 and NGC 7674 show that flie putative jets may 
nevertheless feel an influence from the NLR.) Of the intrinsic 
causes mentioned in the introduction, we cannot constrain 
the role of black-hole spin or magnetic field strength in the 
accretion disc. The expanding plasmon model of Pedlar et al. 
(1985) predicts that the radio components expand into an ISM 
of density 1 cm~^ at speeds initially of > lO'' km s~^ and would 
reach a size of 30 pc after 2000 yr by which time the expansion 



speed would have been reduced to 5000 kms At this speed 
and at the distance of NGC 5506 over the 6.25 yr span of 
our observations, the component diameters would increase 
by 0.3 mas, or, in NGC 7674, over 14.58 yr by 0.1 mas. An 
increase by this amount is, unfortunately, much smaUer than 
could be measured with our data. 

4. Spectral indices: We asked in the introduction whether 
Seyfert nuclei show a single compact synchrotron self- 
absorbed core component at base of a collimated jet, like 
scaled-down radio galaxies. We found flat- or absorbed- 
spectrum components in NGC 5506 and NGC 2110 with 
brightness temperatures of > 3 x 10^ K. The components 
in NGC 5506 have no clear axis of symmetry and the flat 
spectrum could be caused by either free-free absorption or 
synchrotron self-absorption. The morphology of NGC 2110 
shows colUmated jets on VLA scales and so is the only object 
in our sample with a flat-spectrum VLBI core and large-scale 
jets. Thus, the classical conically-expandingBlandford-Konigl 
jet structure (Blandford & Konigl 1979, see also Falcke & 
Biermann 1995) is sometimes but not often seen in our Seyfert 
galaxies. 

5. Statistics of selected radio properties: To improve the re- 
liability of any conclusions we enlarged our sample by drawing 
Seyferts from the literature that were observed with VLBI and 
found that flat- or absorbed- spectrum cores were present in 13 
of 16 cases. The most strongly absorbed components tend to 
be in objects that also show other evidence for absorption, such 
as water maser emission or large X-ray absorption columns. In 
those cases, according to the Seyfert unification schemes, the 
flat-spectrum components would be viewed through an edge- 
on obscuring disk or torus, the ionized inner edge of which 
could be causing free-free absorption. If so, the situation in 
those objects is difl'erent from the Blandford- Konigl jet model. 
However the case for free-free absorption is unambiguous in 
only one galaxy (NGC 3079) for which the spectral index is 
> -h2.5. The other 12 Seyferts have -0.3 < a < +2.5 and 
so spectral turnovers could be produced by either synchrotron 
self-absorption or free-free absorption. 

Where extended structures are present in the enlarged sam- 
ple (14 Seyferts), the position angle diff'erences between pc- 
scale and kpc-scale radio emission were found to be uniformly 
distributed between 0° and 90°. Such bends could be due to 
changes in the jet ejection axis or due to pressure gradients in 
the ISM. No correlation was found between the axis of pc-scale 
radio structure and the rotation axis of the host galaxy. 
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29.7 mas 9.7 mas 
37.1 mas -10.5 mas 




18 4 X 114 




VLBA 1997, 6 cm, Feb 5, 1997 

BO 
Bl 


41.6 ±4.9 
24.4 ± 3.0 


29.6 ± 3.3 
4.3 ±0.7 


0.0 0.0 
28.9 mas 9.1 mas 


0.30 


5.7x3.7 




V 1 . i>/\ lyy 1 , o.D cm, r ell Lyy i 

BO 


27.0 ± 3.4 


21.5 ±2.6 


0.0 0.0 


0.46 


6.5 X 3.0 




FVN IXr^ Nnvll 1Q0Q 

BO 
Bl 
B2 


32.9 ±4.1 
25.6 ± 3.3 
6.1 ± 1.4 


31.0 ±3.5 

22.1 ± 2.6 
5.4 ±1.0 


0.0 0.0 
26.2 mas 9.3 mas 
37.9 mas -15.2 mas 


0.43 


26.2x14.1 


0.06 
-0.73 
-0.30 


MEHLEM, 18 cm, NoY 11, 1999 


113.2 ± 103 


86.6 ± 6.7 


0.0 0.0 


2.41 


259.1x150.9 


-0.33 


FVN (\ rm Fph 28 2111)0 
V 1^, u CHI, f eu ^o, Auuu 

BO 

(tapered) 
Bl 

B2 

(tapered) 


35.9 + 4.1 
32.0 + 3.7 
11.0+ 1.4 
10.8 ± 1.4 
1.9 ±0.5 
4.4 ± 0.7 


36.1 ±4.2 
6.8 ± 1.2 
1.3 ±0.4 


14 13 14.87720 -03 12 27.6460 
14 13 14.87916 -03 12 27.6365 
14 13 14.87982 -03 12 27.6583 


0.20 


8.8x6.8 




MERLIN, 6 cm, Feb 28, 2000 

B0+B1+B2 
(tapered) 


85.8 ± 6.4 
78.8 ± 5.9 


66 8 + 3 8 


0.0 0.0 


0.43 


78 6x67 4 




VLBA 2000, 18 cm. May 31, 2000 

BO 
Bl 
B2 


49.0 + 5.8 
35.8 ±4.4 
4.6 ±1.1 


46.0 + 5.1 
24.7 ± 2.9 


0.0 0.0 
29.7 mas 10.1 mas 


0.45 


18.6x11.3 




VLBA 2000, 6 cm. May 31, 2000 

BO 
Bl 
B2 


48.1 ± 5.2 
15.3 ± 2.7 
2.3 ± 0.8 


42.2 ±4.6 

6.6 ±1.1 

2.7 ±0.7 


0.0 0.0 

29.1 mas 9.6 mas 

45.2 mas -8.5 mas 


0.42 


7.4x4.3 










NGC 7674 








ii>viN, lo Cin, Apr u, i709 

SE 

NW 


31.4 ±3.9 
4.1 ± 1.1 


31.4 ±3.9 
5.1 ±0.97 


0.0 0.0 
-453.1 mas 230.8 mas 


0.46 


39.3 X 24.2 




irVN 1S(<in Nnvll 10O0 

SE 
NW 


37.6 + 4.5 
11.0 + 2.1 


23.9+2.8 
5.6 ± 1.0 


23 27 56.71054 +08 46 44.1310 
23 27 56.67984 +08 46 44.3641 


0.40 


32.5x9.8 


-1.78 


MERLIN 18 cm Nov 11 1999 

SE 
NW 


111.5 ±8.6 
45.6 ±4.0 


74.9 ±4.6 
29.2 ± 2.3 


23 27 56.71730 +08 46 44.0104 
23 27 56.68608 +08 46 44.2195 


0.82 


276x110 


-1.45 
-1.58 


EVN, 6 cm, Feb 28, 2000 

SE 

(tapered) 


6.2+ 1.2 
5.2 ± 1.0 


5.0 ±0.7 


23 27 56.71201 +08 46 44.1368 


0.24 


8.2x4.6 




MERLIN, 6 cm, Feb 28, 2000 

SE 

(tapered) 
NW 


18.0 + 2.3 
187 + 2.4 
3.8 ± 0.7 


16.2 + 1.1 
3.9 ± 0.5 


0.0 0.0 
-452.5 mas 229.5 mas 


0.31 


82.2x36.8 





(tapered) 6.9 ± 1.1 

Table 2. Results from the observations. "Tapered" 6 cm flux densities were measured using images tapered and restored with 
the 18 cm beam size. The MERLIN 18 cm and 6 cm positions differ by 35 mas because the coordinates used at the correlator 
were changed in an unknown way and the resulting offset could not be removed accurately. We have precessed the MERLIN 
B1950 coordinates to J2000 using NRAO's scheduling software SCHED. Positions are absolute for observations that were phase- 
referenced and the phase-referencing worked; such positions are given in h, m, s (for RA) and °, ', " (for Dec). Otherwise the 
positions are given in mas relative to the strongest component. The flux density error estimates represent 1 cr. 
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Int. flux density 


Peak flux density 


Peak position (J2000) 


rms 


Beam size 




(mjy) 


(mjy/beam) 


RA Dec 


(mJy/beam) 


(mas) 


(S OCV") 



EVN,18ciii,N0Yll, 1999 

Core 



MERLIN, 18 cm, Nov 11, 1999 

(tapered) 



Core 


29.4 


+ 


2.9 


North 


11.4 


± 


2.1 


South 


9.1 


± 


1.9 



EVN, 6 cm, Feb 28, 2000 

Core 
(tapered) 



MEJILIN, 6 cm, Feb 28, 2000 

Core 



20.3 + 2.7 
17.6 ±2.4 



13.9 ±1.7 


05 52 11 


.37378 


31.4 + 2.5 


05 52 1 1 


.38263 


12.8 ±2.2 


05 52 11 


.39579 


10.8 ± 2.0 


05 52 11 


,37443 


20.0 ± 2.5 




0.0 


38.9 ± 2.3 




0.0 



VLBA, 18 cm, Aug 6, 1998 

(natural weighting) 



NW 


37.4 ±3.8 


21.6 ±2.5 


08 04 05.85528 


SE 


11.8 ±1.9 


10.9 ±1.8 


08 04 05.85690 


(uniform weighting) NW 


36.0 ± 3.8 


12.8 ± 1.6 


08 04 05.85514 


SE 


12.1 ± 1.9 


8.5 ± 1.5 


08 04 05.85701 



VLBA, 6 cm, Aug 6, 1998 

(tapered) 

NW 

SE 

(full resolution) 



9.1 ± 1.1 
5.0 ± 1.1 



6.3 ±1.2 
4.5 ± 1.0 



3 04 05.85529 +05 06 49.8552 
S 04 05.85694 +05 06 49.8304 



NWl 


1.0 ±0.2 


1.1 ±0.1 


08 04 05.85566 


+05 06 49.8570 


NW2 


2.1 ±0.4 


1.7 ±0.2 


08 04 05.85538 


+05 06 49.8554 


NW3I' 


1.1 ±0.2 


0.9 ±0.1 


08 04 05.85511 


+05 06 49.8553 


NW4'' 


1.2 ±0.3 


1.0 ±0.1 


08 04 05.85500 


+05 06 49.8557 


SE 


3.9 ±0.6 


2.6 ±0.3 


08 04 05.85696 


+05 06 49.8298 


EVN, 18 cm, Nov 11, 1999 










NW 


36.6 ± 4.3 


19.8 ± 2.3 


08 04 05.85520 


+05 06 49.8545 


SE 


13.0 ±2.3 


8.0 ±1.1 


08 04 05.85666 


+05 06 49.8301 


MERLIN, 18 cm, NoY 11, 1999 












121.3 ±9.0 


99.3 ±5.5 


08 04 05.86024 


+05 06 49.8985 


EVN, 6 cm, Feb 28, 2000 












9.8 ± 2.1 


7.7 ± 1.4 


0.0 


0.0 


MERLIN, 6 cm, Feb 28, 2000 












33.1 ± 2.5 


24.1 ± 1.4 


0.0 


0.0 



12.2x8.6 -1.26 
-0.78 



240.7x136.5 -1.01 



(tapered) 



Table 3. (continued) Positions and flux densities measured by fitting Gaussians to the image, instead of fitting a quadratic 
surface and integrating over the component region. 
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Name 


Veron name 


Class 


Lin. Res. 
(pc) 


Instrument 


Ref. 


3C 287.1 




SI 


5.43 


VLBA 5 GHz 


Fomalont et al. (2000) 


3C 390.3 




S1.5 


1.09 


global VLBI 5 GHz 


Alefetal. (1996) 


111 Zw 2 




S1.2 


0.260 


VLBA 15/43 GHz 


Brunthaler et al. (2000) 


Ark 564 




_ 


0.598 


global VLBI 5 GHz 


Vir Lai (2001) 


IC 5063 


PKS 2048-57 


Slh 


3.30 


LBA 2.3 GHz + LBA HI 


Oosterloo et al. (2000) 


MCG 8-11-11 




S1.5 


0.635 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 1 




S2 


6.18 


EVN 1.7 GHz 


Kukula etal. (1999) 


Mrk 1 




82 


0.550 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 3 




Slh 


5.24 


EVN 1.7 GHz 


Kukula et al. (1999) 


Mrk 78 




S2 


1.75 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 231 




Sl.O 


0.343 


VLBA 15 GHz 


Ulvestad et al. (1999b) 


Mrk 231 




Sl.O 


4.09 


VLBA 1.4/2.3/4.8/8.4/15/22 GHz 


Ulvestad et al. (1999a) 


Mrk 273 




S2 


7.32 


VLBA 1.4 GHz 


Carilli & Taylor (2000) 


Mrk 348 




Slh 


0.131 


VLBA 15 GHz 


Ulvestad et al. (1999b) 


Mrk 348 




Slh 


6.41 


EVN 1.4 GHz 


Neff &deBruyn(1983) 


Mrk 348 




Slh 


1.49 


VLBA 1.7/5.0/8/15/22 GHz 


Peck et al. (2003) 


Mrk 477 




Slh 


2.39 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 530 


NGC 7603 


S1.5 


1.45 


global VLBI 5 GHz 


Vir Lai (2001) 


Mrk 766 




S1.5 


0.413 


global VLBI 5 GHz 


Vu- Lai (2001) 


Mrk 926 




S1.5 


0.917 


VLBA 8.4 GHz 


Mundell et al. (2000) 


Mrk 1210 




Slh 


0.392 


EVN 1.6/5 GHz VLBA 1.6/5 GHz 


This paper 


Mrk 1218 




S1.8 


0.582 


global VLBI 5 GHz 


Vir Lai (2001) 


NGC 1052 




S3h 


0.015 


VLBA 5/8.4/22/43 GHz 


Kadler et al. (2002) 


NGC 1052 




S3h 


0.048 


VLBI 1.4-43.2 


Vermeulen et al. (2003) 


NGC 1052 




S3h 


0.038 


VLBA 2.3/8.4/15.4 GHZ 


Kameno et al. (2001) 


NGC 1068 




Slh 


0.221 


VLBA 1.7,5,15 GHz 


Roy etal. (1998) 


NGC 1068 




Slh 


2.94 


EVN 1.4 GHz 


Ulvestad et al. (1987) 


NGC 1167 




S3 


0.895 


global VLBI 5 GHz 


Giovannini et al. (2001) 


NGC 1167 




S3 


0.639 


global VLBI 1.6 GHz 


Giovannini et al. (1990) 


NGC 1275 




S1.5 


0.416 


VLBA 2.3/5/8.4/15.4/22/43 GHz 


Walker et al. (2000) 


NGC 21 10 




Sli 


0.151 


VLBA 8.4 GHz 


Mundell et al. (2000) 


NGC 2110 




Sli 


0.966 


EVN 1.6/5 GHz 


This paper 


NGC 2273 




_ 


0.342 


global VLBI 5 GHz 


Vir Lai (2001) 


NGC 2639 




S3 


0.233 


global VLBI 5 GHz 


Vir Lai (2001) 


NGC 2639 




S3 


0.162 


VLBA 1.6/5/15 GHz 


Wilson etal. (1998) 


NGC 3079 




S2 


0.044 


VLBA 5/8/22 GHz 


Trotter etal. (1998) 


NGC 3079 




S2 


0.054 


global VLBI, 5.0 GHz 


Irwin &Seaquist (1988) 


NGC 3079 




S2 


0.042 


global VLBI 1.4/8.4/15/22 


Sawada-Satoh et al. (2000) 


NGC 3147 




S2 


0.219 


VLBA 1.6/2.3/5/8.4 GHz 


Ulvestad & Ho (2001b) 


NGC 3227 




S1.5 


3.74 


MERLIN 1.6/5 GHz 


Mundell etal. (1995) 


NGC 4151 




S1.5 


0.116 


VLBA 1.6,5 GHz 


Ulvestad et al. (1998) 


NGC 4151 




S1.5 


1.29 


EVN 1.7 GHz 


Harrison et al. (1986) 


NGC 4151 




S1.5 


4.82 


MERLIN 5 GHz 


Pedlar etal. (1993) 


NGC 4168 




S1.9 


0.236 


VLBA 5 GHz 


Nagar et al. (2002a) 


NGC 4203 




S3b 


0.084 


VLBA 1.6/2.3/5/8.4 GHz 


Ulvestad & Ho (2001b) 


NGC 4258 




S2 


0.029 


VLBA 22 GHz 


Hermstein et al. (1997) 


NGC 4258 




S2 


0.232 


VLBA 1.4/1.6 GHz 


Cecil et al. (2000) 


NGC 4395 




SI. 8 


0.221 


VLBA 1.4 GHz 


Wrobel et al. (2001) 


NGC 4565 




S1.9 


0.207 


VLBA 5 GHz 


Falcke et al. (2000) 


NGC 4579 




S3b 


0.118 


VLBA 1.6/2.3/5/8.4 GHz 


Ulvestad & Ho (2001b) 


NGC 4579 




S3b 


0.246 


VLBA 5 GHz 


Falcke et al. (2000) 


NGC 5252 




S2 


0.671 


VLBA 8.4 GHz 


Mundell et al. (2000) 


NGC 5506 




Sli 


0.815 


EVN 1.6/5 GHz 


This paper 


NGC 5548 




S1.5 


0.366 


VLBA 8.4 GHz 


Wrobel (2000) 


NGC 5793 




S2 


0.135 


VLBA 1.7,8.4,15,22 GHz 


Hagiwara et al. (2001) 


NGC 5793 




S2 


0.316 


VLBA+Y27 1.6/5 GHz 


Hagiwara et al. (2000) 



Table 4. 
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Name 


Veron name 


Class 


Lin. Res. 


Instrument 


Ref. 








(pc) 






NGC 5929 




S3 


0.290 


global VLBI 5 GHz 


VirLal (2001) 


NGC 5929 




S3 


6.44 


MERLIN 0.4/1.6/5 


Suetal. (1996) 


NGC 7212 




Slh 


1.37 


global VLBI 5 GHz 


VirLal (2001) 


NGC 7469 




S1.5 


0.844 


global VLBI 5 GHz 


Vii Lai (2001) 


NGC 7674 




Slh 


2.80 


VLBA+Y27+Arecibo 1.4 GHz 


Momjian et al. (2003) 


NGC 7674 




Slh 


2.58 


EVN 1.6/5 GHz 


This paper 


NGC 7674 




Slh 


2.18 


global VLBI 5 GHz 


Vir Lai (2001) 


NGC 7682 




S2 


0.996 


global VLBI 5 GHz 


VirLal (2001) 


T0109-383 


NGC 424 


Slh 


0.226 


VLBA 8.4 GHz 


Mundell et al. (2000) 



Table 4. All Seyferts observed with VLBI. Column 1 lists the names we used, column 2 the names used by Veron-Cetty & Veron (2001), 
column 3 the Seyfert types as described by Veron-Cetty & Veron (2001) , column 4 the linear resolution in pc, column 5 the instrument and 
frequency used, and column 6 the reference to the publication. 



Source 


I'kil/lnv. 


Comp. 


1 X 

< 


Rcfs. 


Mrk231 


* 


N 


-0.86 


Ulvestadct al. (1999a) 






C 


+0.88 


Ulvestad et al. (1999a) 






s 


-1.54 


Ulvestad et al. (1999a) 


Mrk348 


* 




+0.93 


Barvainis & Lonsdale (1998) 


Mrk 463E 




L 


-1.06 


Norris et al., unpublished 






R 


-0.62 


Norris et al., unpublished 






1 


-0.67 


Norris et al., unpublished 






2 


-0.77 


Norris et al., unpublished 






3 


-1.00 


Norris et al., unpublished 


Mrk 1210 




NW 


-1.26 


This paper 






SE 


-0.78 


This paper 


NGC 1068 


* 


NE 


-1.5 


Roy et al. (1998) 






C 


> 2.2 


Roy et al. (1998) 






SI 


> 0.6 


Roy etal. (1998) 






S2 


< -1.9 


Royetal. (1998) 


NGC 2110 




c 


+0.12 


This paper 








1 1 TO 

+ 1. /5 


Wilson et al. (lyys) 


NGC 3079 


* 


E 


+0.77 


Middelberg et al. in prep. 






A 


> +4.54 


Middelberg et al. in prep. 






B 


> +4.36 


Middelberg et al. in prep. 


NGC 3147 


* 




+0.20 


Ulvestad & Ho (2001b) 


NGC 4151 


* 


E 


-0.5 


Ulvestad et al. (1998) 






D 


< -0.2 


Ulvestad et al. (1998) 






F 


< -1.3 


Ulvestad ct al. (1998) 


NGC 4203 


* 




+0.38 


Ulvestad & Ho (2001b) 


NGC 4579 


* 




+0.20 


Ulvestad & Ho (2001b) 


NGC 5506 


* 


BO 


+0.06 


This paper 






Bl 


-0.73 


This paper 






B2 


-0.30 


This paper 


NGC 5793 


* 


CIC 


-0.70 


Hagiwara et al. (2000) 






CINE 


-1.00 


Hagiwara et al. (2000) 






C2C 


-0.46 


Hagiwara et al. (2000) 






C2W 


-0.72 


Hagiwara et al. (2000) 






C2E 


+0.13 


Hagiwara et al. (2000) 


NGC 7469 


* 


2 


-0.34 


Norris et al, unpublished 






3 


+0.09 


Norris et al., unpublished 






5 


-0.12 


Norris et al., unpublished 


NGC 7674 






-1.78 


This paper 



Table 5. Seyferts from Table 4 that have dual-frequency (18 cm and 6 cm) VLBI observations with matching beam sizes, from which the 
spectral index is derived. Asterisks denote objects with at least one component with a flat or inverted spectrum. 
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Sample 




Na>-03 


Motal 


% 


P 


This paper 


our obs + lit search 


17 


38 


45 




Unger et al. (1987) 


X-ray Seyferts 


1 


20 


5 


0.3% 


Rush et al. (1996) 


CfA Seyferts 


3 


46 


6.5 


< 0.1 % 


Rush etal. (1996) 


12 jum Seyferts 


4 


40 


10 


< 0.1 % 


Morganti et al. (1999) 


Distance limited 


13 


60 


22 


1.8% 


Ho & Ulvestad (2001) 


Palomar Seyferts 


21 


45 


47 


50% 


Ulvestad & Wilson (1984) 


Distance limited 


3 


40 


7.5 


<0.1% 



Table 6. Fraction of flat- and inverted-spectrum components in various Seyfert samples. The column 'No-0.3' gives the number of components 
with specttal indices a > -0.3 (i.e., components with flat or inverted spectta). The column Wtotai' gives the total number of components in 
each sample. The column % gives (AfQ,>_o.3/Motai) x 100. The column P gives the level of significance from a comparison of the fraction of flat- 
and inverted-spectrum components in the VLBI sample of this paper (limits treated as detections), to that of each comparison sample at arcsec 
resolution, using the difference-of-two-proportions test with Yates correction for continuity (e.g. Glanz 1992). 









Level of Significance 




Sample 






KS 


Logrank 


Gehan 


Peto-Peto 


Unger et al. (1987) 


X-ray Seyferts 


0.3% 


0.7% 


0.3% 


8.3% 


8.3% 


Rush et al. (1996) 


CfA Seyferts 


0.1% 


0.1% 


0.0% 


0.5% 


0.5% 


Rush ct al. (1996) 


12 fim Seyferts 


0.1% 


0.6% 


0.0% 


2.0% 


2.0% 


Morganti etal. (1999) 


Distance Limited 


0.1% 


0.3% 


0.0% 


1.2% 


1.2% 


Ho & Ulvestad (2001) 


Palomar Seyferts 


1.0% 


27% 


44% 


75% 


77% 


Ulvestad & Wilson (1984) 


Distance Limited 


0.1% 


0.3% 


0.0% 


0.7% 


0.7% 



Table 7. Comparison of the spectral index distribution of the VLBI Seyfert sample to various other Seyfert samples. The level of significance 
gives the probability that the spectral index distribution of the comparison sample and the VLBI seyfert sample in this paper were drawn from 
the same parent population. 
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Name 


Res. 


0(pc) 


(kpc) 


A0 


Array 


Refs. 


© (opt.) 




(pc) 


(deg) 


(deg) 


(deg) 






deg 


3C 390.3 


1.09 
1210 


323 


323 





global VLBI 5 GHz 
VLA 5 GHz 


Alefetal. (1996) 
Alefetal. (1996) 




Mrk78 


1.75 
504 


40 


275 


55 


global VLBI 5 GHz 
VLA 5 GHz 


Vir Lai (2001) 
VirLal (2001) 




Mrk231 


0.343 

3270 


93 


170 


77 


VLBA 15 GHz 

VLA 1.5/8.4/15 GHz 


Ulvestad et al. (1999b) 
Ulvestad et al. (1999a) 


100 


Mrk 1218 


0.582 
294 


330 


317 


13 


global VLBI 5 GHz 
VLA 5 GHz 


Vir Lai (2001) 
Vir Lai (2001) 


130* 


NGC 1052 


0.015 
162 


65 


275 


30 


VLBA 5/8.4/22/43 GHz 
MERLIN 1.4 GHz 


Kadleretal. (2002) 
Kadleretal.(2003) 




NGC 1068 


0.221 
217 


11 


40 


29 


VLBA 1.7,5,15 GHz 
VLA 5 GHz 


Roy et al. (1998) 

Ho & Ulvestad (2001) 


160 


NGC 1167 


0.895 


140 


313 


7 


global VLBI 5 GHz 


Giovannini et al. (2001) 






131 




66 


74 


VLA 5 GHz 


Sanghera et al. (1995) 




NGC 1275 


0.416 

382 


353 


150 


23 


VLBA 2.3/5/8.4/15.4/22/43 GHz 
VLA 1.4 GHz 


Walker et al. (2000) 

Pedlar etal. (1990) 




NGC 3079 


0.044 
310 


125 


60 


65 


VLBA 5/8/22 GHz 
VLA 1.4/4.9 


Trotter etal. (1998) 
Duric et al. (1983) 


75 


NGC 4151 


0.116 
104 


257 


80 


3 


VLBA 1.6,5 GHz 
VLA 5 GHz 


Ulvestad et al. (1998) 
Ho & Ulvestad (2001) 


140 


NGC 4258 


0.029 
125 


3 


355 


8 


VLBA 22 GHz 
VLA 1.5 GHz 


Herrnstein et al. (1997) 
Cecil et al. (2000) 


60 


NGC 7212 


1.37 
609 


257 


25 


52 


global VLBI 5 GHz 
VLA 5 GHz 


Vir Lai (2001) 
VirLal (2001) 


133« 


NGC 7469 


0.844 

193 


28 


353 


35 


global VLBI 5 GHz 
VLA 5 GHz 


Vir Lai (2001) 
Vir Lai (2001) 


35 


NGC 7674 


2.80 
178 


274 


296 


22 


VLBA+Y27+Arecibo 1.4 GHz 
MERLIN 1.7 


Momjian et al. (2003) 
Momjian et al. (2003) 


64" 


nizw2 


0.260 
3460 


303 


240 


63 


VLBA 15/43 GHz 
VLA 1.5 GHz 


Brunthaler et al. (2000) 
Brunthaler, priv. comm. 




Ark 564 


0.598 


282 






global VLBI 5 GHz 


Vir Lai (2001) 


28 


IC 5063 


3.30 
242 


66 


115 


49 


LBA 2.3 GHz + LBA HI 
ATCA 8.3/1 .4GHz 


Oosterloo et al. (2000) 
Morganti et al. (1998) 




NGC 2110 


0.151 
240 


350 


186 


16 


VLBA 8.4 GHz 
VLA 1.4/5 GHz 


Mundell et al. (2000) 
Ulvestad & Wilson (1983) 




NGC 3147 


0.219 


317 






VLBA 1.6/2.3/5/8.4 GHz 


Ulvestad & Ho (2001b) 


65 


NGC 5506 


0.815 
359 


73 


177 


76 


EVN 1.6/5 GHz 
VLA 4.9 GHz 


This paper 
Colbert etal. (1996) 


1 


NGC 7682 


0.996 


125 






global VLBI 5 GHz 


Vir Lai (2001) 


68* 



Table 8. Subset of sources from Table 4 for which both pc and kpc-scale structures have been observed. The misalignment angle, A0, is 
calculated as 0(pc) - ©(kpc), and has then been reflected onto the range of 0° to 90°. For the seven objects at the end of the table, either the 

PA. information of the kpc scale structure was not completely reliable or was unavailable. The position angles of the galaxy minor axes in the 
last colurrm were drawn from de Vaucouleurs et al. (1991), except for those marked a taken from Schmitt & Kinney (2000) and those marked 
b measured by us from Hubble Space Telescope archival images (see §4.3). 



